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ABSTRACT 

We discuss the stability of galactic disks in which the energy of interstellar clouds is 
gained in encounters with expanding supernova remnants and lost in inelastic colli¬ 
sions. Energy gain and loss processes introduce a phase difference between the pressure 
and density perturbations, making disks unstable on small scales for several recipes 
of star formation. This is in contrast to the standard stability analysis in which small 
scale perturbations are stabilized by pressure. In the limit of small scales the dispersion 
relation for the growth rate reduces to that of thermal instabilities in a fluid without 
gravity. If instabilities lead to star formation, then our results imply a secondary mode 
of star formation which operates on small scales and feeds on the existence of a primary 
mode on intermediate scales. This may be interpreted as a positive feedback. Further, 
the standard stability criterion on intermediate scales is significantly modified. 
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1 INTRODUCTION 

Instability and star formation in galactic disks are intimately 
linked. The observational support for that is the significant 
decline of star formation activity in regions with relatively 
low gas mass surface density, as expected from various sta¬ 
bility analyses (e.g. Martin & Kennicut 2001). Theoretically, 
the collapse of cold gas provides ripe conditions for the ap¬ 
pearance of molecular clouds and, subsequently, star forma¬ 
tion. (e.g. Spitzer 1968, Quirk 1972). 

On large scales (of the order of the size of the disk) 
perturbations are stabilized by rotation. Disks are also sta¬ 
ble on small scales if the pressure, P, and density, p, are 
related by P oc p 7 , where 7 is a constant. On intermediate 
scales, neither rotation nor pressure can stop gravity from 
amplifying the perturbations if the surface density is above 
a critical value. The first detailed stability calculations for 
single fluid disks have been done in two classic papers by 
Toomre (1964) and Goldreich & Lynden-Bell (1965). Simi¬ 
lar type of analysis has been generalized to the more realistic 
case of disks containing stellar and gaseous components (Jog 
& Solomon 1984, Rafikov 2001, Griv, Gedalin & Yuan 2002). 
Viscous disks have also been considered and found to be al¬ 
ways unstable, with small scale perturbations growing at a 
low rate that is set by the viscosity (Gammie 1996). These 


studies have assumed that the perturbation in the pressure 
is proportional that in the density. 

Here we consider the instability of the interstellar cloud 
component of the disk, when it is subject to local energy 
gain and loss processes. Clouds gain energy in encounters 
with expanding supernova (SN) remnants, and lose energy 
in inelastic collisions among each other (McKee & Ostriker 
1977, hereafter M077). Both of these processes depend, in 
general, on the local surface density and velocity dispersion 
of the clouds. The star formation rate is enhanced above the 
global average in a region with a positive density fluctuation, 
resulting in a local increase of the velocity dispersion. A per¬ 
turbation in the velocity dispersion changes the energy loss 
rate and may also affect the star formation rate. We per¬ 
form detailed stability analysis under the assumption that 
the cloud component be treated as a hydrodynamical fluid 
over sufficiently large length and temporal scales. The treat¬ 
ment here leads to a significant modification of the standard 
stability criteria. We get instabilities on much smaller scales 
than what is inferred from the standard analysis. In several 
important cases, instabilities extend down to the scale where 
the fluid treatment becomes invalid. 

The paper is organized as follows. The equations are 
derived and expanded in Laplace transforms in §2. Stability 
analysis is presented in §3. This section includes a classifi¬ 
cation of stable and unstable modes by means of Nyquist 
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diagrams, a discussion of a few limiting cases, and a treat¬ 
ment of instabilities with specific form for the energy gain 
and loss rates. The validity of our approach is assessed in in 
55. A summary and discussion of the results are presented 
in § 6 . 


2 THE EQUATIONS 

We write the equations governing the evolution of perturba¬ 
tions in a thin gaseous disk. The gas represents interstellar 
clouds which, for simplicity, are assumed to have isotropic 
velocity dispersion. We assume that the gas has an ideal 
equation of state which after integration over the height of 
the disk yields 

P=(7-1)M£ (!) 


to relate the “two-dimensional” pressure, p, to the mass sur¬ 
face density, p, and the “random” kinetic energy per unit 
mass, £ = Kms/ 2 , where F rms is the three dimensional ve¬ 
locity dispersion of the clouds, and 7 > 1. We will also refer 
to e as internal energy. The adiabatic index 7 relates the 
projected quantities and, therefore, is in general different 
from the physical three dimensional index, T. The relation 
between 7 and T depends on the structure of the disk. For 
a self gravitating disk 7 = 3 — 2/F = 9/5 for T = 5/3 (e.g. 
Garnmie 2001). The difference between the two indices has 
little effect in our applications and so we present numerical 
results only for 7 = 5/3. We write energy conservation in 
the following form, 


de 

dt 


P_dp 
a 2 dt 


+ G — L , 


( 2 ) 


where G(p, e) and L(p, e) are, respectively, the energy gain 
and loss rates per unit mass, assumed to be explicit functions 
of p and s alone. I 11 the meantime we will study some general 
features of the stability of the system without referring to 
specific forms for G and L. 

Let the disk be in a steady state in which the disk is 
axisymmetric, gravity is balanced by rotation, and energy 
gain is balanced by energy loss, i.e. G = L. This steady state 
is characterized by the rotational speed, 12 (r), the gas mass 
surface density, p 0 (r), and the internal energy per unit mass, 
£ 0 (r), as functions of the distance from the center of the disk, 
r. We are interested in the linear response of the system to 
small perturbations from this steady state. For simplicity we 
consider only radial perturbations as the complication in¬ 
troduced by non-axisymmetric perturbations does not hide 
any additional relevant physical effects. Let u, v, p^\ £^\ 
p^\ and <f>^ be, respectively, the perturbations in the ra¬ 
dial velocity, the transverse velocity, the surface density, the 
internal energy, the pressure, and the gravitational poten¬ 
tial field. The first order equations governing the evolution 


of these perturbations are 


P W = (7 - 1)M 0 £ (1) + (7 - l)e 0 M (1) , (3) 

dt £ (1) = $d t p W + C> (1) + C e £ W , (4) 

d t u — 2 Q(r)v + d 7 ■p ( 1 ) // i o( r ) + dr<j>^\z =0 = 0 , (5) 

dtv — 2 B(r)u = 0 , ( 6 ) 

dtp (1) + d r [rp 0 (r)u\ = 0 , (7) 

d r (rd r <t>( v> )/r + d 2 tj>^ = 4nGp 1 <5 D (z) , ( 8 ) 


where p 0 = (7 — l)£ 0 /x 0 , and 




d(G-L) 

dp lf, o’7> 


and C e 


d{G-L) 

d£ “o’ e o ■ 


( 9 ) 


The first two of these equations are obtained by linearization 
of equations (1) and(2), respectively. The equations (5) and 
( 6 ) are the radial and transverse linear versions of the Euler 
equations. Mass conservation is represented by (7) and the 
Poisson equation by ( 8 ). We expand the perturbations in 
Fourier modes, exp(ik ■ r). Considering only modes satisfy¬ 
ing kr 1 , the linear equations lead to, 


i 1 — ttAi + CmM 1 + C e £ 1 (10) 

p 1 + [k 2 — 2nGp 0 k + (7 — l)£ 0 k 2 ]p 1 

+(7- l)M 0 fc 2 £i = 0 , (11) 

where p 1 is defined by p ^ (r, t) = p 1 ( t ) exp(ik ■ r) and sim¬ 
ilarly for £ 1 , and k = 2 f 2 [l + (dln$7/dlnr)/2 ] 1 ^ 2 is the 
epicyclic frequency. In deriving these equations, the rela¬ 
tion 0 r </T^| z= o = i2nGp 1 exp(ik ■ r) has been used (Toomre 
1964). 

The differential equations (10) and (11) are linear with 
constant coefficients and they can be solved by means of 
Laplace transformation (see Appendix A). By taking the 
Laplace transform of (10) and (11) we get 


s£-£i( 0) = (7-1 -Mi(0)) 


+C^M + C e £ , 

( 12 ) 

s 2 M — sp 1 (0) — p 1 (0) + [k 2 — 2nGp 0 k 


+(7 - l)e 0 fc 2 )]A4 + (7 - l)Mo fc2 £ = 0 , 

(13) 


where M (s) and £ (s) are the Laplace transforms of p 1 ( t ) 
and £ 1 ( t ), respectively, and p ± ( 0 ), /q ( 0 ), and £7 ( 0 ) represent 
the initial conditions given at t = 0. Solving for M yields, 


B(s)M(s) = ( s -C £ )[ S ,q(0)+M0)] (14) 

+ (7- l)h o fc 2 [£i(0) - (7- 1 ) 77 -Mi (0)] , (15) 

Mo 

where 

B(s) = s 3 — s 2 C e + suj 0 - C s ojf , (16) 

with 

<*>o(fc) = k 2 — 2irGp 0 k + 7(7 — l)£ 0 fc 2 , (17) 

wi(fc) = wo - ( 7 - l)e 0 fc2 + 7 - lj . (18) 
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The third order polynomial B(s) has at most three distinct 
roots, Sj (j = 1,2,3). Since A4(s) oc 1/B(s) then according 
to the theory of Laplace transorms (see Appendix A) fj , 1 ( t ) 
is a linear combination of exp (sjt). Therefore, an unstable 
(growing) mode corresponds to a root with a positive real 
part. 


ically finding the roots of B(s). Here we discuss a few im¬ 
portant limiting cases that can be studied analytically. Let 
Au ) 2 = u>i — u!q. First consider the limit, |At5 2 -C \u>q\ 
in which the roots differ by small amount from the values 
s = —1 and s = A^J—uS 1 which solve (19) with A ui 2 = 0. In 
this case, 


In the limit of small scales, i.e. large k, the terms pro¬ 
portional to k 2 are dominant in (17) and (18). This means 
that the effects of rotation and self-gravity are negligible. 
In this limit, the characteristic equation B(s) = 0, reduces 
to the dispersion relation derived by Field (1965) for the 
growth rate of thermal instabilities (see Eq. 15 in Field 1965, 
without the thermal conduction term). 


3 STABILITY ANALYSIS 

Our task is to establish the relevant time scales for the evo¬ 
lution of the perturbations and classify the stable and un¬ 
stable modes. The three time scales corresponding to the 
roots Sj of B(s) must be reflected in C e , coo and uq. We will 
find that C e is negative for the specific forms of G and L 
we use below. The dependence on C e can be scaled out by 
working with the variables s = s/{—C e ), u>o = u>o/{—C e ) 
and uq = uJi/(—C e ). If Sj ( j = 1, 2,3) solves 

s 3 + s 2 + suio + aq :s= 0 , (19) 

then Sj = —C e Sj. Therefore, the stability of a mode with a 
wave-number k, depends only on u > 2 (k) and u> 2 (k). The time 
scales, however, may depend nontrivially on C e . 

The roots of B(s) can either be all real (with zero imag¬ 
inary parts) or one real and two conjugate complex roots 
(i.e. having the same real part). The sum of the roots is 
si + S 2 + S 3 = C e <0 and therefore at least one of them 
must have a negative real part, meaning that there is always 
a decaying mode. The question of the existence of unstable 
modes for any combination of uiq and lo 2 can be addressed 
by Nyquest diagrams (Nyquist 1932, see also Appendix B). 
These diagrams provide a simple way to determine the num¬ 
ber of roots lying to right of the imaginary axis (i.e. roots 
with positive real parts) without having to solve the third 
order equation B(s) = 0. The outcome of this analysis is 
as follows (we refer only to roots lying to the right of the 
imaginary axis), 

(i) lOq > to 2 > 0 : no roots 

(ii) uif > ujo > 0 : 2 complex roots 

(iii) u>q > 0 , uif < 0 : 1 root 

(iv) u>q < 0 , u>l < 0 : 1 root 

(v) Uq < 0 , u > 2 > 0 : 2 root (either 2 real or 2 complex 
conjugates) 

Therefore, the disk is stable only for uig > ui 2 > 0. 


3.1 limiting cases 

Nyquest diagrams do not contain any information on the 
instability time scales. Those could be obtained by numer¬ 


Sl 


S2,3 


I 1 + ( LOx/Ce ) 2 

|Oe| l + (a, 0 /C e )2 ’ 

| ±|Ce| + w 0 U>1 ~ U>o /—— 

|a| — CTT^ -2^r v_u;o 


( 20 ) 


Note that these expressions involve C e in a complicated 
way. For lo\ > w 2 > 0 the system is unstable, in ac¬ 
cordance with the Nyquist analysis, with a growth time 
scale of [Re(s 2)] _1 = [(1/2)C e (w 2 — u;o)/(Cf +Wo)] -1 - For 
u>q > ujf > 0 there are no unstable modes. 

The second and more important limit is when cl>o 1 
and ACj 2 = 0(u!q). This holds for small scale perturbations 
since for sufficiently large k the leading term in w 2 is the one 
involving k 2 (see eq. 17). The roots of B(s) are approximated 
by, 

2 2 2 

Sl = -\C e \'^2 and S 2,3 = | Ce \ ——— !f° ± lUJo ■ ( 21 ) 

o > 0 2 uJ 0 

For u ! 2 > uj 2 the root si represent a decaying mode, while 
the conjugate roots S 2 and S 3 correspond to sound waves 
that grow at a rate proportional to |C e | and oscillate with 
frequency of wo/2tt. The origin of the growth is heating oc- 
curing during the compression the waves (c.f. Field 1965). 
In terms of C M and C e the condition, w 2 > u>g, for the ap¬ 
pearance of growing sound waves reads, 


Mo Cm 

£ 0 C e 


+ 7 - 1 < 0 , 


which is similar to the condition in Eq.24 of Field 1965. 
For u>q > uj 2 > 0 all modes are decaying with S 2 and S 3 
representing overstable (oscillating) modes. For uj 2 > 0 but 
u > 2 < 0 the mode si is growing without oscillations and S 2 
and S 3 are overstable modes. The mode described by si cor¬ 
responds to growth by condensation under nearly constant 
pressure conditions as described by Field (1965). It can re¬ 
sult from an enhanced cooling efficiency as the density is 
increased. The condition uif < 0 holds for 


2 _ MqC/j 

£ 0 C e 


< 0 . 


( 22 ) 


This is equivalent to the condition for a condensation mode 
as given by Eq.23 of Field 1965. 


3.2 Numerical solutions 

To better visualize the time dependence of perturbations we 
plot in figure ( 1 ) the functions /q(t) and ^(t) for u > 2 < 0 
(top panel), uj 2 > uj 2 > 0 (middle) and u > 2 > uj 2 (bottom). 
The figure shows only cases with u > 2 > 0 which is always true 
for sufficiently large k. These curves have been obtained by 
numerical integration of the following dimensionless form of 
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( 10 ) and ( 11 ), 


~ u . 

Pi + 


1 — -p—k + k I jl 1 + k^E 1 — 0 

Q o 

K = (7 - l)A'i + i + C e £ 1 


(23) 

(24) 


where the prime symbol denotes differential with respect 
to the variable tn, ko = nGfi 0 /(('y — l)e 0 Qo), k = k/ko 
Qo = k[( 7 - 1 )e 0 ] 1/2 /( 7r GM 0 ), = e 1 /e 0 , /q = /q//z. 0 > 

C M = yU 0 C' M /(e 0 /K), and C e = C s /k. All solutions are for 
C e = —1, k = 5, Qo = 1, and 7 = 5/3. The solutions in the 
top, middle and bottom panels correspond to C M = —1.3, 
0, and 2, respectively. The initial conditions of all solutions 
are p, 1 (0) = e 1 (0) = 1 and /i' = 0. The behavior of the 
solutions is in agreement with Nyquest diagrams and the 
limiting cases. For u)f > u>g > 0 (bottom panel), the solution 
is oscillatory with a growing envelope. These growing sound 
waves arise because of the slight heating during compression. 
The case u>q > u> 2 > 0 is overstable. Note that the derivative 
of e 1 at r = 0 is positive although the dotted curve may 
appear as declining at r = 0. The phase difference between 
the dotted and solid curves, however, becomes more evident 
as r increases. 

For u>i < 0 and u> 2 > 0 the solution is unstable without 
oscillations. This mode is driven by thermal instability of a 
condensation mode (Field 1965). Note the phase difference 
between /q and e 1 . This phase difference is the drive behind 
the instabilities shown in the top and bottom panels. 


4 STABILITY WITH SPECIFIC FORMS OF G 
AND L 



0 1 2 3 4 5 6 

T — t/C 

Figure 1. The time evolution of the perturbations in the surface 
density (solid line) and internal energy (dotted) in arbitrary units. 
The curves have been obtained by direct numerical integration of 
the equations of motion. 


A key parameter is the ratio ^i 0 C IJ ,/£ 0 C s which determines 
the difference between lo\ and u>o (see Eq.18). This ratio 
depends on the specific forms chosen to describe the energy 
gain and loss functions. 

The function L(/x, e) represents energy dissipation per 
units mass, via inelastic cloud collisions. Hereafter we adopt 
the following form (M077, Efstathiou 2002), 

L = rj^fj^E 1 ' 2 . (25) 

The energy gain function, G, is directly proportional 
to the star formation rate (M077). We work with two dis¬ 
tinct choices, the first is a Schmidt-Kennicutt (e.g. Kenni- 
cutt 1989) form and the second assumes that the star forma¬ 
tion rate depends explicitly on the Toomre parameter which 
governs instabilities on intermediate scales. 


4.1 Stability with a Schmidt-Kennicutt star 
formation law 

Here we write the energy gain function as, 

G = Va» n , (26) 


where r/ G is a numerical factor having the appropriate units. 
Equating this G with L given by (25) gives, 

C M = (n - 2 ^c/hT 1 and — ^o~" • ( 27 ) 

z Vg 

Therefore, using (18), 

u> 2 = uig + (7 — l)(2n — 3 — 7 )s 0 k 2 . (28) 

Of particular interest is the small scale behavior. For 
sufficiently large k we write u)g « 7(7 — 1 )e 0 fc 2 . Therefore, 
for n < 3/2, we have ui 2 < 0, while for n > ( 3 + 7 )/ 2 , we have 
Wi > u>g. The disk is unstable in these two cases. But the 
nature of the instability is different: the former is described 
by a single non-oscillatory mode, while the latter has two 
growing oscillatory modes (see fig. 1 and the discussion of 
Nyquist diagrams at the end of §3). 

We plot in figure (2) the real parts of the roots as a 
function of the wave number k for n = 1.5 (top panel) and 
n = 2 (bottom). The dashed line is the instability growth 
rate according to the standard analysis a la Toomre (1964). 
This is given by \J—u>q for Uq < 0, and zero otherwise. A 
value of 7 = 5/3 is used and n, ^0 and £ o have been tuned 
so that the maximum of the dashed curve is consistent with 
the growth rate inferred from the lower curve of Fig.la of 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

A _1 [kpc _1 ] 


Figure 2. The growth rate of the unstable modes as function of 
the scale A -1 = k/2ir. The top and bottom panels show results 
for two values of the power index of the star formation rate, n. 
The solid line is the real part of the root of B(s). The dashed 
line is the growth rate estimated as (—u>q) 1//2 for ojg < 0 and 
zero otherwise. The values of the various parameters were chosen 
such that the minimum of u)g matches the value given in Jog & 
Solomon (1984). 


Jog & Solomon (1984). We use \C e \ = 20 km s _1 kpc -1 , 
corresponding to a time scale of 5 x 10' yr (M077). The solid 
line shows the positive real parts of the roots. According to 
the solid curve the system with n = 1.5 remains unstable 
on small scales, albeit with a longer time scale than the 
instabilities at intermediate scales. For n = 2 the instability 
extends to smaller scales than what would be inferred from 
the standard analysis (dashed line). 


4.2 Q dependent star formation rate 


Assume now that the star formation rate depends on 
the time scale defined by standard stability analysis (e.g. 
Toomre 1964). We work with a gain function of the form 
(e.g. Wang & Silk 93, Elmegreen 1999), 

G = n G KnF(Q) (29) 


where 


Q 


ne 1 ' 2 

nGfj, 


and 


p (1 - Q 2 ) 1/2 

F - Q -• 


(30) 


Note that the meaning of r/ G is different from (26); the same 
symbol is used only for the sake of brevity. Energy balance 
(G=L) at ft = fj, 0 and e = e 0 yields 


F(Q) = Q/K , 


(31) 



Q 


Figure 3. The ratio uj^ /lUq as a function of the Toomre parameter 
for the Wang & Silk star formation rate. The ratio is computed in 
the limit k — > oo. The system is unstable when the ratio exceeds 
unity. 


where 1/K = 


L )fi 2 Q . This relation gives, 


Q =- 


K 2 K 2 / 4 

2 + 2 V 1 + K 2 


1/2 


(32) 


For K 1 and If > 1 we have Q = K 1//2 and Q = 1, 
respectively. The first order variation of (G — L) is 


SG — 8L — rj G t 


„ AF dQ AF dQ 
^ Ml+M dQ^ 1 +P dQ^ £l 

- 2 ?? L /ie 1/2 /ii - l») L fi 2 e _1/2 e 


(33) 

(34) 


Comparing this expression with the definitions of C^ and 
C e (see Eq. 9) and using dQ/d^i = —Q/n, dQ/de = Q/{2e) 
and £ 1,/2 = 2 -kGhQ/k we get, 


Ce 


~Vg k 



d In F \ 

d In Q J 


F, 


-Vg 1 


Mo 

2e„ 



d In E \ 

d In Q J 


F . 


(35) 

(36) 


Since din E/d In Q < 0, the coefficient C e is negative. By 
evaluating the logarithmic derivative in terms of Q we find, 
aC^/eCe = 2Q 2 /{Q 2 - 2 ), and 

~ (7 ~ 1 )e 0 fc2 ^7 — 1 + ■ (37) 

For Q 2 > 2 ( 7 —l)/( 7 +l), we have u> 2 > u>g, yielding unstable 
modes even on small scales. For lower q, the system is stable 
on small scales but, as in the case with n = 2 of the previous 
subsection, unstable modes extend to scales smaller than in 
the standard analysis. The ratio uj 2 /l 0 q as a function of Q is 
plotted in figure (3), in the limit of very large k. A value of 
7 = 5/3 is used in this figure. 


5 VALIDITY OF THE APPROACH 

The treatment of the clouds as a hydrodynamical fluid 
is valid over scales larger than the mean free path, Zf p , 
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for cloud-cloud collisions. In the solar neighborhood Zf p « 
100 pc for the cold neutral cores and about 10 pc if the warm 
envelops are included (M077). Although the warm envelops 
should contribute to the collisions, we take Zf p « 100 pc as 
an upper limit. 

Further, there are two time scale to be considered. 
First, there is the mean time between collisions, t co ii = 
Zfp/Kms ~3x 10 6 yr for V^ms = 10 kms -1 . Second, there 
is the mean time, t re heat, between two successive “heating” 
events of a cloud. This comes about because we assume 
that multiple encounters with with expanding SN remnants 
eventually amount to an increase in random kinetic energy 
rather than bulk motions. According to Cox & Smith (1974) 

Zreheat ~ 10® - 10 7 yr. 

There is also the length scale over which the velocity dis¬ 
persion can equilibrate in a time t = max(t re h ea t,Icon) This 
scale arises from the finite mean free path of the clouds 
(see also Gammie 1996) and is similar to the usual heat 
conduction (e.g. Zel’dovich & Raizer 2002). We estimate 
this scale as (tKmsZff,) 1 ^ 2 . Taking t = f co ii yields a scale 

Of (tcollFrmsZfp) 1 ' ~ Zfp. 

Therefore, for an environment like the solar neighbor¬ 
hood, our approach is valid over scales larger than a hundred 
parsecs. The scale is sufficiently small that our treatment re¬ 
mains interesting. For example, this scale is smaller than the 
scale of unstable modes found by means of standard stability 
analysis (e.g. Toomre 1964, Jog & Solomon 1984). 

Finally, there is a continuous process of cloud destruc¬ 
tion and production. We assume that this process is rapid 
and always produces a fixed spectrum for the cloud size dis¬ 
tribution. 


6 SUMMARY AND DISCUSSION 

We have considered the instability of gaseous disks subject 
to local energy gain and loss processes in the presence of 
star formation. The gas represents interstellar clouds and 
the energy is gained in repeated encounters of the clouds 
with expanding SN remnants and is lost in inelastic cloud- 
cloud collisions. These energy exchange processes introduce 
a phase difference between the density and pressure pertur¬ 
bations. In several interesting situations this phase difference 
causes the pressure to amplify the density perturbation. For 
a star formation rate proportional to /x" the instabilities ex¬ 
tend to much smaller scales than what is inferred from stan¬ 
dard stability analyses a la Toomre (1964) and Goldreich & 
Lynden-Bell (1965). 

The small scale instabilities may be responsible for trig¬ 
gering further star formation. Therefore, there may be two 
modes of star formation, a primary mode on intermediate 
scales and a secondary mode operating on small scales at a 
rate that is determined by the coefficient C e (see Eq. 9). This 
can be interpreted as positive feedback in which the interme¬ 
diate scale mode is driving star formation on smaller scales 
through the development of instabilities. It is unclear what 
gain function one should use in (2). A fully self-consistent 
stability analysis should incorporate the energy gain result¬ 


ing from the star formation mode that is induced by the in¬ 
stabilities. This complicates the problem substantially since 
the form of the gain function in this case must be derived 
from the stability analysis self-consistently. Nevertheless, we 
have found that small scale instabilities develop at some 
level for most generic forms of the gain function. We do not 
expect a full self-consistent treatment to change our conclu¬ 
sions significantly. 

We have not included the coupling of gas to the stellar 
component. This could easily be done, but does not affect 
the main conclusions of the present work. Including the stel¬ 
lar component would tend to destabilize small scale modes 
even further (Jog & Solomon 1984). 
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APPENDIX A: A BRIEF REVIEW OF 
LAPLACE TRANSFORMS 

We briefly review some of the relevant properties os Laplace 
Transforms The Laplace transform, /(s), of a function F(t), 
where t > 0, is defined as 

POO 

/(s) = C{F(t)} = / exp(-st) F(t)dt. (Al) 

Jo 
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We will need the Laplace transforms of first and second 
derivatives of a function. Using (Al) these transforms can 
be related to the f(s) by 
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£{F\t)} = sf(s)-F( 0) 

C{F"{t)} = s 2 f(s) — sF(0) — F'(0) , (A2) 

where the prime and double prime denote first and second 
order derivatives, respectively. The Bromwich integral ex¬ 
presses F(t) in terms of f(s) as 

i /»7+ioo 

F(t) = — / exp (st) f(s)ds (A3) 

J 'y—ioo 

where i = \J — 1 and 7 is a real number chosen so that all 
poles of f(s) lie, in the complex plane, to the left of the 
vertical line defining the integration path. Therefore, by the 
residue theorem we have 

F(t) = [residues of exp (st) f(s)] (A4) 

As an example consider f(s) = l/(s — si)(s — S2) which 
has two simple poles at s = si and S2 . The residues 
of exp (st)f(s) at these poles are exp(sif)/(si — s 2 ) and 
— exp(« 2 t)/(si — S2) so that, by (A4), F(t) = [exp(sif) — 
exp(s 2 t)]/(«i — 52). If si = S2, the function has a pole 
of order two at si. The residue in this case is d[(s — 
si ) 2 exp(st)/(s))]/ds evaluated at s = si. Therefore, F(t) = 
t exp(sit). 


APPENDIX B: STABILITY BY MEANS OF 
NYQUIST DIAGRAMS 

Given a polynomial, Pn, of degree, N, Nyquist diagrams 
provide an easy and elegant method for determining the 
number of its roots which lie to the right of the imaginary 
axis, i.e., roots with a positive real part. The method can be 
summarized as follows: a) plot a contour of the value of the 
polynomial in the complex plane when its argument, s, is 
varied from s = +ioo to s = —too along the imaginary axis 
and b) use the contour diagram to compute the increase, 
A ip, in the phase of Pn(s) as s moves along the imaginary 
axis. The number of roots with positive real parts is then 
(Ai/) + Nn)/(2n). 

Figure (Bl) shows Nyquist diagrams for the polynomial 
B(s) given in (16) for three cases as indicated in the figure. 
In the top panel the phase changes from — 7 r /2 to — 37 r /2 
and therefor the number of roots with positive real parts 
is 1. Further, in this case this root has zero imaginary part 
since otherwise its conjugate is also be a root and there 
would be two roots with positive (and equal) real parts. In 
the middle panel the change is — 37 r and no roots lie to the 
right of the imaginary axis. In the bottom panel there are 
two roots with positive real parts. 



Re(B) 


Figure Bl. Nyquist diagrams of the polynomial B(s). 


© 0000 RAS, MNRAS 000, 000-000 




